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Abstract 

We report a study of the suppressed decay B~ [K^ 'k^]£,K^ [and its charge-conjugate mode) 
at BeUe, where [K~^t:~]d indicates that the K~^'k~ pair originates from a neutral D meson. A data 
sample containing 274 million BB pairs recorded at the T(4S') resonance with the Belle detector at 
the KEKB asymmetric e'^e" storage ring is used. This decay mode can be used to extract the CKM 
angle 1^3 using the so-called Atwood-Dunictz-Soni method. The signal for B~ [K~^tt~]dK~ has 
2.7(7 statistical significance, and we set a limit on the ratio of B decay amplitudes rs < 0.28 at the 
90% confidence level. We observe a signal with 5.8a statistical significance in the related mode, 
B- [K+7T-]Dn~. 

PACS numbers: 
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INTRODUCTION 



The extraction of 03, an angle in the Kobayashi-Maskawa triangle[l], is a challenging 
measurement even with modern high luminosity B factories. Several methods for measuring 
03 use the interference between D^K~ and B^ D^K", which occurs when 

and decay to common final states[2]. In this paper, we analyze the suppressed decay 
B^ —>■ [K^n^juK' and its charge conjugate mode, where [K'^ir']]:, indicates that the K~^'ir~ 
pair originates from a neutral D meson. In this case, the color-allowed B decay followed 
by the doubly Cabbibo-suppressed D decay interferes with the color-suppressed B decay 
followed by the Cabbibo-allowed D decay(Fig.l). This decay mode can be used to extract 
03 using the so-called Atwood-Dunietz-Soni method(ADS method) [3]. 




FIG. 1: B [K+TT-jDK- decays. 



ADS METHOD 



Here we define the amplitudes for B decays and D decays as follows: 



Ab = A{B- 
Ab = AiB- 



The branching fractions for B^ 
states / are given as follows: 



T{B- 
T{B- 



{f]DK- 



where 



AiB+^D'K-'), Ao^A{D'^f)=A{D'^f), 
A{B+^D'K+), AD = A{D'^f)=A{D'^f). 

[fjoK^ decays with and Z)° decays to common final 



[r| + rl + 2rBrD cos(-03 + 6)]\Ab\^\Ad\ 
[r| + rl + 2rBrD cos(03 + 6)]\Ab\^\Ad\^ 



6 = 6b + Sd 





Ab 




Ad 


tb = 


Ab 


td = 


Ad 



and 5b and 5d are the strong phase differences between the two B and D decays, respectively. 
The modulus of the amplitude, lA^p can be measured using a flavor specific decay mode. 
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If we use a D decay mode in which Ad and arc known, the above 2 equations have 3 
unknowns(03, tb, S). However, using two final states /i and /2, there are 4 equations and 4 
unknowns(03, tb, Si, 62), which can be solved for 03. Using multiple decay modes for D —>■ fi, 
the value of ^3, and other unknowns, can be extracted from a fit. The suppressed decay 
B~ — i> [K~^71~]dK~ is an especially useful mode for the ADS method. The two interfering 
amplitudes in this decay mode are comparable, and large CP violating asymmetries can be 
expected. This decay mode is thus sensitive to the value of ^3. 

ANALYSIS 

In this paper, we report an analysis of the suppressed decay — > [K^7r^]r)K^. We 
also analyzed the suppressed decay B"^ — > [K^Tr^J^jTr^. In addition, the allowed decays 
B^ — * [K^7!'^]r)K'^ and B^ [K^TT^]D'n'~ are used as control samples to reduce systematic 
uncertainties. The same selection criteria for the suppressed decay modes arc applied to the 
control samples whenever possible. Throughout this report, charge conjugate states are 
implied except where explicitly mentioned and we denote the analyzed decay modes as 
follows. 

Suppressed decay B~ — > [K'^7r~]Dh~ : B~ — > D^sh^ 
Allowed decay B~ [ir-7r+]c/i" : B' Dfh-{h = K, tt) 

The results are based on a data sample containing 274 million BB pairs, collected with the 
Belle detector at KEKB asymmetric energy e~^e~ collider operating at the T{4S) resonance. 
The Belle detector is a large-solid-angle magnetic spectrometer that consists of a sihcon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of timc-of-fiight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals located 
inside a super-conducting solenoid coil that provides a 1.5 T magnetic field. An iron fiux- 
return located outside of the coil is instrumented to detect Kl mesons and to identify muons 
(KLM). The detector is described in detail elsewhere [4]. Two diff'crent inner detector 
configurations were used. For the first sample of 152 million BB pairs, a 2.0 cm radius 
beampipe and a 3-layer silicon vertex detector were used; for the latter 122 million BB pairs, 
a 1.5 cm radius beampipe, a 4- layer silicon detector and a small-cell inner drift chamber were 
used [5]. 

Event selection 

D mesons are reconstructed by combining two oppositely charged tracks. These charged 
tracks are required to have a point of closest approach to the beam line within ±5 mm of 
the interaction point in the direction perpendicular to the beam axis(dr) and ±5 cm in the 
direction parallel to the beam axis{dz). A K/n likelihood ratio P{K/ti) = Ck/{jCk + ^n) is 
formed for each track, where Ck and are kaon and pion likelihoods. We used the particle 
identification requirement P{K/'k) > 0.4 and P{K/7r) < 0.7 for kaons and pions from 
D Ktt decays, respectively. D candidates are required to have an invariant mass within 
±2.5(7 of the nominal mass: 1.850 GeV/c^ < M{K7r) < 1.879 GeV/c^. To improve the 
momentum determination, tracks from the D candidate are refitted according to the nominal 
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D° mass hypothesis and the reconstructed vertex position (a mass-and-vertex-constrained 
fit). 

B mesons are reconstructed by combining D candidates with primary charged hadron 
candidates. For the charged tracks, we require P{K/ti) > 0.6 for the kaon in B~ — > 
DK^ and P{K/Tr) < 0.2 for the pion in Dn~ . The signal is identified by two 

kinematic variables, the energy difference AE = + Ex-(n-) ~ -E-bcam and the beam- 
energy-constrained mass Mbc = y^-E'beam ~ (Pd +pW-(7r-))^) where E^ is the energy of the 
D candidate, Ek-(tt-) is the energy of the i^~(7r~) and i?beam is the beam energy, in the cm 
frame, pr, and PK-{n-) are the momenta of the D and K~{7r~) in the cm frame. We define 
the signal region as 5.27 GeV/c^ < Mbc < 5.29 GeV/c^ and -0.04 GeV < AE < 0.04 GeV. 
In the case of multiple candidates per event, we choose the best candidate on the basis of 
a determined from the difference between the measured and nominal values of Md and 
Mbc. 



qq continuum suppression 

To suppressed the large background from the two-jet like e'^e" — > qq{q — u, d, s, c) 
continuum processes, variables that characterize the event topology are used. We construct 
a Fisher discriminant of Fox- Wolfram moments called the Super- Fox- Wolfram(5'FW^) [6] [7], 
where the Fisher coefficients are optimized by maximizing the separation between BB events 
and continuum events. Furthermore, cos 9b, the angle in the cm system between the B flight 
direction with respect to the beam axis is used as another variable to distinguish BB events 
from continuum events. These two independent variables, SEW and cos 9b are combined 
to form a likelihood ratio(Li?), 

LR — /^sig/(-^sig ~l~ -^cont) 

p _ fSFW y /'COS 0s 

''-'sig(cont) — ''-'sig{cont) ^ ''-'sig(cont) ' 

where £sig and £cont are likelihoods defined from SEW and cos 9b distributions for signal 
and continuum backgrounds, respectively. We optimized the LR requirement by maximizing 
a figure of merit, S/^/S + N, where S and denote the expected number of signal and 
background in the signal region. For B^ DcsK~ {ti~) we require LR > 0.85(> 0.75), which 
retains 44.8%(57.6%) signal events and removes 96.2%(93.2%) of the continuum background. 



Peaking backgrounds 

For B~ — > DcsK~, one can have a contribution from B~ — > £)°7r~, D'^ — > K'^K~ , which 
has the same final state and can peak under the signal. In order to reject these events, 
we veto events that satisfy 1.843 GcV/c^ < M{KK) < 1.894 GeV/c^. The allowed decay 
B~ Dfh~ can also be a peaking background for the suppressed decay modes due to Kn 
misidentification. Therefore, we veto events for which the invariant mass of the Ktt pair 
is inside the D mass cut window when the mass assignments are exchanged. Furthermore, 
three-body charmless decays B~ K~^K~7t~ and B~ K^ti~ti~ can peak inside the 
signal region for B^ D^sK^ and B^ D^sT^^ , respectively. These peaking backgrounds 
are estimated from the AE distributions of events in a D mass sideband, defined as 1.808 
GeV/c^ < MiK-K) < 1.836 GeV/c^ and 1.893 GeV/c^ < M(ir7r) < 1.922 GeV/c^, which 
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are shown in Fig. 2. For B~ — > Dcs7r~ , the peaking background estimated by fitting the 
plot is consistent with zero. Since the Standard Model prediction for the K^7T~7r~ 
branching fraction is smaller than 10~^^[8], this background contribution is ignored. On the 
other hand, for B~ — > DcsK~, the estimated peaking background is 3.1 ± 2.9 events inside 
the AE signal region after scaling to the D mass signal region. As a check, we naively 
estimate the expected background from the measured B~ K^K~tt~ mode. According 
to [9], the B- K+K-TT- yield is 94 ± 23 events with an efficiency of 13.8% (78.7 fb"^). 
Using this result, the estimated background is 

257.1 fb"^ areao bSdk ^„ 

94 X -j- X X — ~ 2.9 events, 

78.7 fb Sire&Daiitz gSkk^ 

where we assumed that the B~ — > K^K~ti~ yield is uniformly distributed over the Dalitz 
plot, and area/j/area/ja/jf^ is the ratio between the D mass cut area and the Dalitz plot 
area, and Q^DK/Q^KKn{= 17.5/13.8) is the ratio of the B~ — > DcsK~ efficiency to the 
B~ — > K'^K^'K" efficiency. This naive estimate is consistent with the estimate from the 
mass sideband. Therefore, we subtract 3.1 ± 2.9 events from the observed B" — > DcsK~ 
yield. 

After applying all the cuts, the signal efficiencies are 17.5% and 24.6% for B~ DcsK~ 
and B~ —>■ Dcs7t~ , respectively. The signal yields are extracted by fitting the AE distribu- 
tions. 



Fitting the AE distributions 

Backgrounds from decays such as B~ Dp~ and B~ D*tt~ are distributed in the 
negative AE region and make a small contribution to the signal region. The shape of 
this BB background is modeled as a smoothed histogram from generic Monte Carlo (MC) 
samples. The continuum background populates the entire AE region. The shape of the 
continuum background is modeled as a linear function. The slope is determined from the 
AE distribution of the Mf,c sideband data (5.20 GeV/c^ < M^^ < 5.26 GeV/c^). 

The AE fitting function is the sum of two Gaussians for the signal, the linear function 
for the continuum, and the smoothed histogram for the BB background distribution. 

In the fit to the AE distribution of B~ — > Df7r~ , the free parameters are the position, 
width and area of the signal peak, and the normalizations of continuum and BB back- 
grounds. The ratio of the two Gaussians of the signal is fixed from the signal MC. For the 
B~ — > DfK~ fit, the position and width of the signal peak are fixed from the B~ — > Df7r~ 
fit results. To fit the feed-across from -D/7r~, we use a Gaussian shape where the left and 
right sides of the peak have different widths since the shift caused by wrong mass assignment 
makes the shape asymmetric. The shape parameters of this function are fixed at values de- 
termined by the fit to the B~ — > Dfii^ distribution using a kaon mass hypothesis for the 
prompt pion. The areas of signal and feed-across from Dtt^, and the normalizations of con- 
tinuum and BB backgrounds are fioated in the fit. For B~ DcsK~ and B~ — > Dcs7^~ , the 
signal and BB background shapes are modeled using the fit results of the B~ D jK" and 
B~ —>■ DfT:~ modes, respectively. The area of the feed-across from Dcs7t~ is estimated as 
the measured yield of B~ — > DcsT^~ multiplied by the n to K misidentification probability. 
However, the areas of the signal and the normalizations of continuum and BB backgrounds 
are fioated. The fit results are shown in Fig. 3. The numbers of events for B~ — > Dcshr 
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and Dfh , and the statistical significances of the B — > Dcsh signals are given in Table 

I. The statistical significance is defined as 1^—2 ln(£o / '^max) ? where /^max is the maximum 
likelihood in the AE fit and Cq is the likelihood when the signal yield is constrained to be 
zero. The uncertainty in the peaking background contribution is taken into account in the 
statistical significance calculation. The statistical significance of the B" — > DcsT^~ signal is 
over 5.0(7. 

TABLE I: Signal yields and efficiency. For the B~ — ^ DcsK~ signal yield, the peaking background 
contribution has been subtracted. 

Mode Product branching Efficiency Signal Yield Statistical 

fraction from PDG (%) significance 

B- ^ DcsK- - 

B- DcsTT- (6.9 ± 0.7) X 10-^ 

B-^DfK~ (1.4 ± 0.2) X 10-5 

B- — DjTT- (1.9 ± 0.1) X 10"^ 



17.5 lb 0.2 U.rtji 2.7 

24.6 ± 0.2 30.7111 5.8 
17.5 ± 0.3 535.0l}|^ 

2 1.7 ± 0.2 10178t|n I 



RESULTS 

Branching fraction of suppressed decay modes 

The branching fractions for i?^ — > D^sh^ih = K^ti) are determined as 

B{B- ^ D,,h-) = B{B- ^ Dfh-) X 

J-^Dfh 

where Nd^^h and Nofh are the number of B~ — > Dcsh~ signal events and B~ — > Dfh~ 
signal events. The product branching fractions for B~ Dfh~ , calculated from the world 
averages for the branching fractions [10], are given in Table 1. Using these, the branching 
fractions for the suppressed decays B~ — > Dcsh~ are found to be 

B{B- [K+7r-]DK-) = {3.9tll{stat) ± 0.2{sys) ± O.Q{PDG)) x 10"^ 
B{B- [K+7r-]D7r-) = {5.7t\-l{stat) ± 0.3{sys) ± 0.3{PDG)) x 10"^ 

Most of the systematic uncertainties from the detection efficiencies and the particle iden- 
tification cancel when taking the ratios, since the kinematics of the B~ — > Dcsh~ and 
Dfh~ processes are similar. The systematic errors arc due to the uncertainty in the 
yield extraction and the efficiency difference between B~ — > Dcsh~ and B~ — > Dfh~. The 
uncertainties in the signal shapes and the qq background shapes are determined by varying 
the shape of the fitting function by ilcr. The uncertainties in the BB background shapes 
are determined by fitting the AE distribution in the region -0.07 GcV < AE < 0.20 GeV 
ignoring the BB background contributions. The uncertainties in the efficiency differences 
are determined by the signal MC. The total systematic errors are obtained as the quadratic 
sum of those uncertainties. The results are shown in Table II. 
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The uncertainties in the branching fractions arc statistics-dominated. For the B — > 
DcsK~ branching fraction, we set an upper hmit at the 90% confidence level as 

B(B- D^,K-) < 7.6 X 10-^(90% C.L.), 

where we took the likelihood function as a single gaussian with width given by the quadratic 
sum of the statistical and systematic errors, and the area is normalized in the physical region 
of positive branching fraction. 

TABLE II: Systematic uncertainties for B~ — ^ DK~ and B~ — ^ Dir" . 







Systematic 


error (%) 


Source 




- DfK- 


Dcsvr" -D/vr" 


BB background shape 


±2.1 


±1.0 


±4.6 ±1.6 


qq background shape 


±3.6 


±0.4 


±1.9 ±0.1 


Signal shape 


±0.6 


±0.4 


±1.4 ±0.2 


Feed-across shape 


±1.4 


±1.0 




Efficiency difference 




±1.5 


±1.3 


PDG NormaUzation 


±14.3 


±5.3 


Total 


±4.9 ± 14.3(PBG) 


±5.5 ± 5.3(pi3G) 



Ratio of branching fractions Rdh 



We define the ratio 



^""^ = B{B- Dfh-) + B{B^ Dfh^) - ^'^^ 



Dfh 

The ratios Rdh are determined as follows 



Rdk = {2.8tli{stat) ±0.1{sys)) x 10"^ 
Rd^ = {S.Otolistat) ±0.2{sys)) x IQ-^ 



and 



Rdk < 4.7 x 10-^(90% C.L.). 
The ratio Rdk is related to 03 by 

Rdk = + + 2rBr£) cos 03 cos 5, 

where 



td 



A{D^ K-7r+) 



= 0.060 ± 0.003. 
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Using the above result, we obtain a limit on tb- The least restrictive limit is obtained 
allowing ilo" variation on rn [10] and assuming maximal interference (03 = 0°, 5 = 180° or 
03 = 180°, 5 = 0°) and is found to be 

Tb < 0.28. 



CP asymmetry 

We search for partial rate asymmetries Aoh in — > Dcsh^ decay, fitting the and 
B~ yields separately for each mode, where Aoh is determined as 

^ B{B- ^ D,sh-) - B{B+ D,sh+) _ ^ ^ 
- B{B-^D,sh-)+B{B+^D^h+) ^^-^'^^ 

The peaking background for B~ — > DcsK~ is subtracted assuming no CP asymmetry. The 
fit results are shown in Fig.4 and Table III. We find 

Adk = OAQf'oilistat) ±OM{sys), 
Ad^ = 0.12toi7{stat)±0m{sys), 

where the systematic uncertainty is from the intrinsic detector charge asymmetry, the B~^ 
and B~ yield extraction, and the asymmetry in particle identification efficiency of prompt 
kaons. The intrinsic detector charge asymmetry is determined from the B"^ — > Dfi:^ sam- 
ples. The systematic uncertainty from yield extraction is determined by varying the fitting 
parameters by ilo". The systematic uncertainty due to particle identification efficiency of 
prompt kaons is explained in [12]. The total systematic errors are combined as the quadratic 
sum of those uncertainties (Table IV) . The measured partial rate asymmetries Adh are con- 
sistent with zero. 



TABLE III: Signal yields and partial rate asymmetries. 



Mode 


N{B-) 


N{B+) 


Adh 


B ^ D,sK 

B ^ DcTT 


11 2+6-1 
-'■■'-•^-5.4 

-1 7 9+6.5 

-■^'•^-5.8 


13.61^1 


0.49l[j:|| ± 0.06 
0.12l°;^? ± 0.06 



TABLE IV: Source of systematic uncertainties for the asymmetry calculation. 





Systematic 


error (%) 


Source 


Adk 


Ad-k 


Yield extraction 


4.8 


4.9 


Intrinsic detector charge asym 


2.5 


2.5 


PID efficiency of prompt kaons 


1.0 




Total 


5.5 


5.5 
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SUMMARY 



Using 274 million BB pairs collected with the Belle detector, we report studies of the 
suppressed decay B" Dcsh~{h = K,7r). We observe B" Dcs'k' for the first time, 
with a significance of 5.8o". The size of the signal is consistent with expectation based on 
measured branching fractions [10]. The significance for B~ DcsK~ is 2.7a and we set 
an upper limit on the ratio of B decay amplitudes tb- This result is consistent with the 
measurement of tb in the decay B~ — > DK~, D — > Kstt^t^' [13]. 
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FIG. 2: A.E distributions for events in the mass sideband for B DcsK (left) and B 
-DcsTT" (right). The signal shapes are modeled using the results of the B~ — > Dfh~{h = K,'k) 
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FIG. 3: AE fit results for B' DcsK-{top-Mt), B' DcsTt" (top-right), B' D/K" (bottom- 
left), and B~ Djtt" (bottom-right). The charge conjugate modes are included for these plots. 
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FIG. 4: AE fit results for B' Dcs^" (top-left), B+ Dcs^"^ (top-right), B' DcsTt" (bottom- 
left), and B'^ — DcsTt"^ (bottom-right) . 
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